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Motivated by recent experiments on low-dimensional quantum magnets in applied magnetic fields, 
we present a theoretical analysis of their properties based on the nonlinear a model. The spin 
stiffness and a 1/N expansion are used to map the regimes of spin-gap behavior, predominantly 
linear magnetization, and spin saturation. A two-parameter renormalization-group study gives the 
characteristic properties over the entire parameter range. The model is relevant to many systems 
exhibiting Haldane physics, and is applied here to the two-chain spin ladder compound CuHpCl. 



PACS numbers: 75.10.Jm, 75.30. Cr, 75.40. Cx 

The importance of low-dimensional spin systems in re- 
vealing fundamentally new quantum mechanical proper- 
ties has been recognized since Haldane's conjecture 0] 
concerning the effects of quantum fluctuations in integer- 
and half-integer-spin antiferromagnetic (AF) chains. We 
begin with the experimental observation that the mag- 
netization curves of some such materials, thought to be 
prototypical of the extreme quantum limit, are in fact re- 
markably classical. As examples, we refer here to the Hal- 
dane (S = 1) chain NENP 0, and primarily to the two- 
chain 5=1 ladder CuHpCl ||J|, which show regimes 
of linear magnetization whose gradient is the (classical) 
Neel susceptibility. 

With the goal of understanding such behavior, we con- 
sider the quantum AF system in an external magnetic 
field using the nonlinear a model (NLsM). This widely- 
applied treatment is in fact semi-classical, being truly 
valid only in the limit of large spin S, but has in the past 
formed the basis for many fundamental deductions con- 
cerning the quantum limit. We will demonstrate here its 
applicability for effectively integer-spin quantum systems 
in appreciable magnetic fields, and provide justification 
of this result in terms of the suppression of quantum fluc- 
tuations by the field. For the purposes of developing the 
present theoretical analysis, we will concentrate on the 
best-characterized sample in recent literature, Cu2(l,4- 
diazacycloheptane^CLi (CuHpCl) HQ], 

The Hamiltonian for the ladder system in a magnetic 
field b = .g/^sB may be written as 



H = 
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m,i '^m,i+l 



(i) 



where J is the intrachain exchange interaction and J' the 
interchain, or ladder "rung" , interaction. The deriva- 
tion of the NLsM in the presence of a magnetic field 
is presented in Ref. M. For A^-site chains with pe- 
riodic boundary conditions along x, in the geometry 
shown in Fig. 1, the resulting model has J x = J, 
J y = ij', b = (0,0,6), and S m ,jv J+i = S mji . There 
are two key points. First, the spin is written in terms of 
slowly- varying, orthogonal, staggered and uniform com- 



the lattice constant, and the fluctuations 1 about the 
staggered configuration must be integrated out subject 
to the orthogonality constraint n • 1 = 0. Second, the 
full Euclidean action in space and inverse temperature, 
Se = Sb + Jq drH, contains in addition a Berry-phase 
term S B = iA J drdx[-al-(n A h)] + AmS(Pi + P 2 ). 
Pi = P2 = j- J drdx(n A ii)-d x n separate to give simply 
the Pontryagin index on each chain when d y n = 0, as is 
the case in a ladder of only two chains. The last term 
in Sb is thus i(4irPi)2S, demonstrating that the system 
will have integer-spin characteristics for any value of S, 
and the topological term may be ignored ||,|7j . 

The action for the quasi-one-dimcnsional ladder sys- 
tem, in 1+1 Euclidean dimensions denoted by /i, is then 



Se = ^ / P M n) 2 - (b 2 - (n • b) 2 ) + 2^b•n A h], 

*9 Jr.x 



(2) 



: (2/NyS)\fj/Jx is the bare coupling constant, 
It 



where g 

J = Jx + ^Jy, and the integral over r is to upper limit 
Lt = cJ3, with c = 2Sa\/ J X J /% (h — > 1) the effective 
spin-wave velocity. We have left explicit the number N y 
of chains in the ladder; for the compound under consid- 
eration J y > J x and N y — 2, giving an effectively rigid 
rung coupling. The form of the NLsM in an external 
field given by the second term in Eq. (||) has been de- 
rived previously S, and its implications considered in 
the low-field limit. In what follows we will examine its 
effects for arbitrary fields. 

To gain initial insight into the effect of the magnetic 
field, we consider the spin stiffness of the ladder system 
||. Taking the staggered spin configuration to be subject 
to a twist 9 in the plane normal to the applied field (Fig. 
1), we calculate the free energy F(b,6) to one-loop or- 
der in the important out-of-plane spin fluctuations, and 

deduce the spin stiffness from p s = \cL QJ? 

z att 6=0 



E 



LL T ^ k 2 + (6/c) 2 



(3) 



ponents as S r 



S[(-l 



air, 



where a is where p° = c/2g is the classical value and the sum 



1 



provides both quantum and thermal (through the finite 
"length") corrections to first order. 
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FIG. 1. Representation of two-chain ladder system, show- 
ing directions of staggered moment, twist and applied field. 

We restrict the analysis to the low-temperature, or 
"quantum" case Lt > L. Evaluating the summation 
between spatial limits tt/L and tt / 'a, and introducing the 
"magnetic length" L m as n/L m — b/c, 

■{a/Lf + {a/L m f 
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(1 + (a/L m ) 2 ) 

L in Eq.(||) may be regarded as the correlation length 
£, beyond which segments of the ladder behave inde- 
pendently, and can be computed from p s = 0. In 
the weak-field limit (L m — > oo) one recovers the result 



Co 



e a7r5 , where a = y/J x /J. The 



general solution is £ (B) = £o/ \Jl — (L* yl /L m ) 2 , where 
L* m = aVe 4 */9 - 1 gives the critical field B* at which 
the correlation length diverges. For fields B < B* , the 
finite correlation length may be written as £,(B) = ae a ^ s , 
where S = S[l - ,g/47rln(l - (L* m /L m ) 2 )] is a growing 
value of the effective spin. For B > B* , the field enforces 
a quasi-long-ranged correlation throughout the system 
||, and it is most convenient to write the spin stiffness 
as p s = p°[l — g/47rln(l + (L m /a) 2 )], which recovers the 
bare value as B — > oo. Finally, the divergence of the cor- 
relation length at B* corresponds to the closing of the gap 
A to spin excitations according to A a y/l — (B/B*) 2 . 
This situation is summarized in Fig. 2. 

However, if the B derivative is taken of the free en- 
ergy in order to compute the magnetization, the result- 
ing terms do not yield the required zero result in the spin 
gap regime. This indicates the breakdown of the approx- 
imation implicit in deriving F(b, 9) that the field be large 
on the scale of in-plane fluctuation energies \<j)\. Thus the 
spin stiffness analysis, while qualitatively revealing of the 
behavior of the system, is not reliable at low fields. 

The situation in the weak-field regime may be ad- 
dressed by a 1/N expansion Jl(|, which is expected to 
be appropriate in describing spin-gap phases. The stag- 
gered spin n is taken to exist in an TV-dimensional spin 
space, in which only the component n z is selected by the 
field. The relevant parts of the O(N) action are 



U(n 2 



1) 
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where b denotes b/c and the constraint that n have unit 
magnitude is made explicit with the Lagrange multiplier 
iX, whose saddle-point value is given by 



1 



k 2 + iX + b 2 ' 



(0) 



in which the B-field term is found to appear only at 
0(1/ N) . iX functions as a mass, or cutoff term in momen- 
tum integrations, and is thus an upper lengthscale for co- 
operative processes in the system, or simply a correlation 
length (inverse excitation gap). Writing the saddle-point 
solution as iX — c 2 ir 2 /£(_B) 2 and carrying out the sum- 
mation at low T gives an expression analogous to p s = 
emerging from Eq.(|J). 




FIG. 2. Schematic depiction of behavior of correlation 
length or spin gap, and of spin stiffness, with applied field. 

At weak fields (£ < L m ) one finds £ ~ ae 2rr / N s(l + 
(a 2 /NL 2 n )e 47r/g , a 0(1/ N) correction to a result differ- 
ing from the previous one by a power of 1/N in the ex- 
ponent. The result Jg (^Ip - ) = ~~W ensures both that 
the corresponding magnetization contribution from the k 
summation terms in F is identically zero to 0(1/ N), and 
that the constant term b 2 + c 2 n 2 /\ 2 yields 26(1 - 1/N). 
It is clear that the behavior required of a gapped system 
is returned in the weak-field regime only on making the 
well-recognized identification (motivated by comparison 
with renormalization-group results 0]) N -> N-2, and 
by returning to the physical situation N — 3. Then the 
magnetization is indeed zero, and the saddle-point solu- 
tion for £,(B) becomes precisely that deduced from the 
spin stiffness analysis, with the same critical field B* . 
With this replacement caveat we thus obtain from these 
two approaches a consistent picture of both weak- and 
strong-field regimes. 

To quantify the regimes of validity of the foregoing 
analyses, we consider next a renormalization-group (RG) 
approach to the NLsM in an applied field. Taking the 
model in the form 



S, — / dxKc^n) 2 - ■/,-(] 
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and transforming to variables <f> and ^fgo z = n z repre- 
senting respectively the in- and out-of-plane fluctuations, 
the latter in a form suitable for perturbative expansion 
in the coupling constant g, we obtain 



2g K > 2 y ' 



A)a z + 0(g). (8) 
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A denotes the in-plane terms (d^(f>) 2 + 2ib(p, which be- 
cause <f>(r, x) is assumed to vary slowly can be taken to 
be a small constant (no fast Fourier modes) in the mo- 
mentum shell 77r/a < |k| < ir/a (7 — > 1). Performing 
the integral and expanding in A, the form of Eq.(||) is 
recovered with coefficients g{a') and b(a') 2 given by par- 
tial traces. Evaluation of these and differentiation with 
respect to the flow parameter I — ln(a'/a), leads to the 
coupled RG equations 



dl 
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1 + /3 2 ), (9) 



in which (3 — a'b(a'). These new RG equations possess a 
variety of interesting limiting cases, whose detailed study 
we defer to a future publication Jl^ |. 

For the present purposes, we concentrate on the fixed 
points to obtain a qualitative picture of the RG flow 
diagram, and on the consequences for the magnetiza- 
tion, (i) Seeking a fixed point by weak-field expan- 
sion around (3* = 0, we find dg/dl ~ g 2 /2n and thus 
d\nj3 2 /dl = 2 — d\ng/dl, which may be solved to yield 
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The fixed point (g*,/3*) = (00, 0) is clearly stable if the 
flow is stopped at I* = 2ir/go. The system will flow 
to this strong-coupling regime if the starting value j3~o is 
sufficiently small. The lengthscale £* = ae 1 ' at which the 
flow stops may be compared with the spatial and thermal 
dimensions L, Lt of the system to calculate directly the 
effects of finite size and temperature ]l3| . (ii) At strong 
fields (/3 — > 00), dg/dl — 0, or j = go , indicating that 
the coupling is not renormalized, and dln/3 2 = 2dl, from 
which it follows that b(l) — bo, i.e. neither is the field. 

Numerical solution of Eqs.@ leads to the flow diagram 
in Fig. 3. The regime (ii) of strong initial field may be 
termed the weak-coupling situation, where g and b are 
weakly renormalized to finite values. Here the perturba- 
tion theoretic approach is consistent and the weak cou- 
pling corresponds to deconfincment of the excitations on 
the lengthscale L m set by the field. The regime (i) is the 
strong-coupling condition, of confinement of (gapped) ex- 
citations, where the assumption of small g in the deriva- 
tion proves to be inconsistent. However, in this regime 
one may deduce the critical lengthscale £* governing the 
behavior of the system and, as we will show below, that 
the magnetization is zero (inset Fig. 3). The critical 
starting field separating the two regions is &* = 0.46. 
The properties of the two regimes may be illustrated by 
considering the correlation length in each. £ is a phys- 
ical quantity and does not change under the RG flow, 
meaning that d^/da' = 0, whence 
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(i) At small p 2 , — <9£/£ ~ 2irdg/g, which has solution 
£0 = ^e 271 ^ 1 / 30-1 / 9 '. Under the reasonable assumption 
that £ — > a, the bare lattice constant, as the system flows 
to the strong-coupling limit (1/go ~ > 0), we recover the 
expected result £0 = ae 27T / go (see above) for the finite 
physical correlation length, (ii) For large (3 2 , —2d^/^ ~ 
d(3 2 /f3 2 leads to £0 = (,£/ a (bo invariant) under the RG 
flow, and we see that for any assumed finite £, the bare 
correlation length £ is the system size L or Lt- 




FIG. 3. RG flow diagram for g and b. Strong- and 
weak-coupling regimes are separated by separatrix s. 

We return to the experimental motivation for the above 
analysis, and compute the magnetization M — OF/ dB of 
the model over the full field range, 



M = -~g^ B N x N y (b/4J) + M d + M 



(12) 



M contains contributions linear in B from the quadratic 
term (H), Md from the dynamical term in the trace (c/. 
@) due to fluctuations of n out of the plane perpen- 
dicular to B, and Mf from in-plane fluctuation terms 
4>. The last giv e a sawtooth form leading to magneti- 
zation steps [pL 3| , a finite-size effect which will not be 
considered further here. While the linear term is al- 
ways present, we have shown that below a threshold field 
B* , where the system has a spin gap, it is cancelled by 
the corresponding correlation-length term. Above an up- 
per threshold B C 2, the magnetization will saturate at the 
value M s = g/j,BSN s (N s — N x N y ), and this effect is not 
contained in the (large-S") model as applied. Systematic 
inclusion of an additional total spin constraint is possible 
Jl3[ | , but to compress the analysis we will here apply M s 
as a simple cutoff. The "dynamical" contribution is given 
consistently to lowest order in the small parameter c/b, 
and in limit of low T, by the constant Md — ^N x gfig. 
Next-order corrections involving the spin excitations have 
the form M' d - Bln(B/J). 

Specializing to the two-chain ladder material CuHpCl, 
the exchange constants deduced from the magnetization 
and susceptibility ||§ are J' = 13. 2K and J = 2.4K, 
whence J = 13. 3T and J x = 3.6T. Taking the simplest 
case of constant Md, and the lower critical field B c \ for 
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onset of magnetization where M{B C \) — 0, we obtain 
B c i = J/g = 6.6T. The saturation field B c2 is given from 
M(B c2 ) = M s as B c2 = AS/gJ = 13.3T. These values 
are in remarkably good agreement with a linear extrap- 
olation of the magnetization data at lowest temperature 
in Rcf. If , where B cl = 6.8T and B c2 = 13. 7T. Such an 
extrapolation appears closer to the data than the predic- 
tions |ll| of a repulsive boson model. We note further 
that the gradient of the linear magnetization is precisely 
the classical Neel susceptibility \af — (iJ^b) 2 /(AJa 2 ) 
per unit volume, a result in itself remarkable for an AF 
in the quantum limit. The computed magnetization is 
shown in Fig. 4, where the dashed line indicates the va- 
lidity limit of the calculation. 
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FIG. 4. Computed magnetization of ladder system CuHpCl. 

The predicted magnetization shows a surprising degree 
of accuracy for a model strictly valid only in the limit 
of large S. We note that the NLsM has recently been 
used with considerable success also for small systems with 
small S |Q . Its accuracy in the current context may be 
ascribed in part to the suppression of quantum fluctu- 
ation effects by the field, as shown above, such that in 
the high- field regime, which covers much of that beyond 
B c \, the behavior of a system with no topological term 
is rather classical. Weak interladder interactions caus- 
ing the real material to display 3d order at intermediate 
fields will also enhance classical behavior. 

Further perspective on the above results is given by 
computing the critical field B* from the spin stiffness 
and 1/N treatments. We deduce from L* m that B* = 
0.32 J/g = 2. IT, a value rather lower than B c \ above, 
and emphasize that in the regime between these values, 
as also seen in the RG analysis, the treatment is not 
reliable. In the RG approach, the magnetization 
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from which we see in the small- B regime, where 
dg/dfto = and d/3/dj3 — e l y/l — gol/^n, that for any 
F(b) analytic in b, M(B ) ~ y/l ~ g Q /2nln(C/a) -> as 



C — > C* . Thus the magnetization vanishes in the strong- 
coupling limit as required. From the numerical solution, 
the field scale 6* gives B* = 1.6T, in reasonable agree- 
ment with the value B* above. We have already shown 
that the regime approaching the physical B c \ lacks a suit- 
able calculation scheme. Physically, the problem is one 
of determining the effects of quantum fluctuations of the 
spin system when the spin gap is weakened by the applied 
field, and the quantum disordered phase thus made less 
robust. We have detailed above the lowest-order predic- 
tions of the NLsM for the closing of the gap, and expect 
these to be appropriate when the gap becomes smaller 
than other energy scales (quantum, thermal and finite- 
size fluctuations) in the system. 

The foregoing analysis is not restricted to CuHpCl, but 
applies also to the S = 1 AF ("Haldane") chain. NENP 
is considered to be a prototypical case for observation 
of the Haldane gap but for the complication of a large 
single-ion anisotropy. The present study predicts again 
the qualitative features of a gapped regime of zero mag- 
netization, followed by approximately linear behaviour 
towards saturation (not achieved), as in experiment [0. 
We will present elsewhere the quantitative aspect of this 
problem. There has been considerable recent interest in 
the possibility of magnetization plateaus in certain sys- 
tems, and we observe that in the current model these may 
be expected, for example in S > 1 chains, when the field 
strength is such that the projected in-plane spin Sn± is 
of integer amplitude, leading to a gapped phase. 

In summary, the nonlinear a model treatment repro- 
duces well the behavior of quantum antiferromagnets in 
an external field. For effectively integer-spin systems, 
meaning those with a trivial topological term, this state- 
ment is valid even in the extreme quantum limit of low 
spin and low dimensionality, because the magnetic field 
acts to suppress quantum fluctuation effects. 
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